REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  inclixling  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed  and  completing  and  reviewing  this  collecUon  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  Information,  Including  suggestions  for  reducing 
this  burden  to  Department  of  Defense.  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arilngtori,  VA  22202- 
4302  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  e  collection  of  Information  if  It  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS.  _ _ _  - 

i.  REPORT  DATE  (DD^MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

1-03-2000  Journal  Article  1999 

A  TitiF  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


r  REPORT  DATE  (DD^MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  ( 

X-03-2000  Journal  Article  1999 

4.  TITLE  AND  SUBTITLE  /  5a.  CONTRACT  NUMB 

Change  in  microscopic  tunneling  on  a  macroscopic  time  scale: 

Current  surge  model  . 

5b.  GRANT  NUMBER 


6.  AUTHOR(S) 

Danhong  Huang,  DA  Cardimona,  and  Anjali  Singh 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force “Rese^xch  Laboratory 
Space  Vehicles  Directorate 
3550  Aberdeen  Ave.  SE 

Kirtland  AFB,  m  87117-5776  _ 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5c.  PROGRAM  ELEMENT  NUMBER 

62601F 


5d.  PROJECT  NUMBER 

4846 


5e.  TASK  NUMBER 

CR 


5f.  WORK  UNIT  NUMBER 

Cl 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  is  Unlimited. 


Z005020T  on 


13.  SUPPLEMENTARY  NOTES  ^  - ' 

Published  in  Journal  of  Applied  Physics,  Volume  87,  Number  5,  pages  2427  -  2430,  1  Mar  2000. 


14.  ABSTRACT 

The  current  Surge  model  is  proposed  to  calculate  the  tunneling  current  in  multiple  quantum  well  structures  in  the 
presence  of  a  slow  time-dependent  electric  field.  The  microscopic  origin  of  a  zero-bias  residual  current  observed 
and  reported  previously  [A.  Singh  and  D.  A.  Cardimona,  Opt.  Eng.  (Bellingham)  38, 1424  (1999)]  is  explored.  The 
mystery  of  the  observation  of  a  microscopic  tunneling  change  on  a  macroscopic  time  scale  is  uncovered,  and  an 
involvement  of  a  very  slow  physical  process  is  shown.  Some  new  features,  such  as  a  current  “ripple,”  a  current 
instability,  a  current  hysteresis,  and  a  current  “arch,”  are  predicted  and  confirmed  experimentally. 


15.  SUBJECT  TERMS 

hysteresis,  quantum  wells,  microscopic. 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Mr.  David  Cardimona 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Unlimited 

5 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(505)  846-5807 

standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


I 


VOLUME  87,  NUMBER  5 


1  MARCH  2000 


^  -  TOURHAL  OF  APPLIED  PHYSICS 

Change  of  microscopic  tunneling  on  a  macroscopic  time  scale: 
Current  surge  model 

Danhong  Huang,®*  Anjali  Singh,  and  D.  A.  Cardimona 
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The  current  surge  model  is  proposed  to  calculate  the  tunneling  current  in  multiple  quantum  weU 
structures  in  the  presence  of  a  slow  time-dependent  electric  field.  The  microscopic  origin  of  a 
zero-bias  residual  current  observed  and  reported  previously  [A.  Singh  and  D.  A.  Cardimona  Opt 
Eng.  (Bellingham)  38,  1424  (1999)]  is  explored.  The  mystery  of  the  observation  of  a  microscopic 
tunneling  change  on  a  macroscopic  time  scale  is  uncovered,  and  an  involvement  pf  a  very  slow 
physical  process  is  shown.  Some  new  features,  such  as  a  cunent  “ripple,”  a  current  instability,  a 
current  hysteresis,  and  a  current  “arch,”  are  predicted  and  confirmed  experimentally  ©  2000 
American  Institute  of  Physics.  [S0021-8979(00)01205-6] 


I.  INTRODUCTION 

In  a  recent  article,*  the  two  authors  found,  experimen¬ 
tally,  a  residual  tunneling  current  in  multiple  quantum  well 
(MQW)  structures  when  a  time-dependent  bias  voltage, 
which  steps  firom  —  5  to  5  V  in  steps  lasting  firom  5  to  60  s, 
sweeps  through  zero.  A  circuit  model,  which  introduced  an 
important  leakage  resistance  in  series  with  a  capacitance, 
was  devised  to  explain  this  phenomena,  and  a  satisfactory 
numerical  simulation  was  obtained.  This  simulation  indicates 
that  a  physical  process  with  a  very  large  time  constant  is 
involved.  However,  the  microscopic  origin  of  this  observa¬ 
tion  remained  a  mystery.  It  is  known  that  the  QW  capaci¬ 
tance  in  MQW  structures  is  on  the  order  of  Cqw~  10"  *°  F  , 
leading  to  a  negligible  charging  energy  e^HCfgf,  when  an 
electron  is  added  to  a  QW.  This  is  in  contrast  to  the  Coulomb 
blockade  effect^  observed  in  quantum  dots  which  results  in  a 
much  smaller  capacitance  compared  to  Cq^.  However,  the 
resistance  of  the  sequential  electron  tunneling  (SET)  in  the 
MQWs  can  reach  as  high  as  R — 10**  fi  at  low  bias  voltages 
and  temperamres.  Consequently,  the  macroscopic  charging/ 
discharging  time  (MCDT)  in  the  MQWs  is  in  the  range  of 
10  s,  which  is  much  greater  than  the  microscopic 
resonant  electron  tunneling  (RET)  time  r,<  1  ps  in  a  double 
QW.  The  question  that  arises  is:  Why  can  we  observe  a 
change  of  the  microscopic  tunneling  on  a  macroscopic 
charging/discharging  time  scale?  In  order  to  resolve  this  is¬ 
sue  and  to  get  a  complete  understanding  of  the  microscopic 
origin  of  the  previously  observed  zero-bias  residual  tunnel¬ 
ing  current  (ZBRTC),  as  well  as  to  provide  in-depth  theoret¬ 
ical  explanations  for  some  recently  observed  phenomena 
which  we  describe  as  a  current  “ripple,”  a  current  instabil¬ 
ity,  a  current  hysteresis,  and  a  current  “arch,”  we  introduce 
a  current  surge  model  (CSM).  Our  proposed  model  success¬ 
fully  reproduces  and  explains  all  of  the  observations. 
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It  is  obvious  that  RET  can  occur  only  when  the  barrier 
between  the  adjacent  QWs  is  thin.  If  the  barrier  is  very  thick, 
the  phase  of  the  wave  function  will  be  lost  as  an  electron 
tunnels  from  one  barrier  to  another.  As  a  result,  only  SET 
exists  for  thick  barriers.  When  the  electric  field  £),(/)  is  fixed 
at  zero,  the  SET  current  /a[f),(/)]  is  also  zero.  However, 
will  be  finite  when  4(r)  is  nonzero  and  varies 
extremely  slowly  with  time.  In  this  case,  electrons  in  the 
MQWs  are  said  to  be  in  a  steady  state,  and  /a[i),(/)]  is 
driven  by  i),(r)  only  through  nonlinear  current-voltage  char¬ 
acteristics.  If  £j(f)  is  switched  suddenly  from  one  value  to 
another  at  a  particular  time,  /«[£),(/)]  will  deviate  from  the 
previous  steady  state  and  sjowly  relax  into  a  new  steady 
state.  We  would  like  to  determine  the  characteristic  time  for 
^a[^*(^)]  to  switch  from  the  old  to  the  new  steady  state.  Of 
course,  within  the  adiabatic  model,  4[5t(t)]  can  follow 
fj,(t)  instantaneously,  switching  from  one  steady  state  to  an¬ 
other  with  no  time  delay.  Therefore,  in  order  to  find  the  true 
switching  time,  we  have  developed  the  CSM  that  moves  well 
beyond  the  adiabatic  limit. 

The  organization  of  this  article  is  as  follows.  In  Sec.  n, 
we  present  our  current  surge  model  beyond  the  adiabatic 
limit  for  deriving  the  dynamical  equation  of  the  nonadiabatic 
field  in  the  system.  Numerical  results  and  discussions  are 
presented  in  Sec.  HI  for  the  total  electron  sequential  tunnel¬ 
ing  currents  under  either  a  step-like  or  a  sinusoidal  time- 
dependent  electric  field.  Experimental  results  and  compari¬ 
sons  are  shown  in  Sec.  IV  for  the  total  electron  sequential 
tunneling  current  under  a  sinusoidal  time-dependent  bias 
voltage.  The  article  is  concluded  in  Sec.  V. 

II.  MODEL  AND  THEORY 

Let  us  start  the  CSM  by  considering  first  the  steady-state 
transport  of  electrons  in  the  MQWs  in  the  presence  of 
If  Shit)  is  time  independent,  the  electrons  undergo  a 
Fowler-Nordheim  tunneling  process'*  through  the  barrier  to 
give  rise  to  a  dc  current.  Here,  the  charge  density  (CD)  in 
each  QW  remains  at  its  equilibrium  value  «2d  •  If  S^ft)  be- 
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comes  time  dependent,  the  electrons  respond  to  it  by  produc¬ 
ing  an  ac  current.  If  at  each  time  depends  on 

instantaneously,  this  tunneling  process  is  said  to  be  adiabatic, 
and  will  not  remember  an5d:hing  that  happened  earlier.  In 
reality,  however,  we  find  that  the  existence  of  current  hyster¬ 
esis  indicates  that  tunneling  electrons  do  remember  what  oc¬ 
curred  previously,  and  thus  the  tunneling  process  must  be 
nonadiabatic  in  nature.  In  fact,  we  find  that  the  tunneling 
current  depends  not  only  on  but  also  on  d£i^{t)ldt.  In 
the  nonadiabatic  tunneling  process,  the  CD  in  each  QW 
changes  with  time  and  fluctuates  around  n2D  when  is 
applied,  which  is  described  well  by  a  nonadiabatic  field 
This  field  can  be  either  positive  or  negative  when  the 
CD  in  the  QWs  is  lower  or  higher  than  /Z2d  ♦  In  the  follow¬ 
ing,  we  will  establish  a  dynamical  equation  for  within 
the  CSM. 

Using  Levine’s  SET  current  formula,^  we  are  led  to 

(1) 

where  is  the  electron  temperature  and  S  is  the  sample 
cross-sectional  area.  The  drift  velocity  5^,(0]  in  Eq.  (1)  is 
related  to  £b(t)  by  a  saturation  model,^  and  the  effective 
three-dimensional  electron  tunneling  density  nQf££i^(t),Tg'] 
is  defined  through 

r+oo 

dET[E,£,(t)] 

JO 

X{fo[E+E,]-ME+Ei  +  eLs£t,(t)-]}, 


where  jidTe)  is  the  chemical  potential  of  electrons  in  each 
QW  and  /o(Z)={l+exp[X-/t,(r,)/fc5r,)F'. 

is  independent  of  5^(0  and  is  denoted  sim¬ 
ply  hy  El  —  /jLc(T^),  In  Eq.  (2),  m*  is  the  effective  mass  of 
electrons,  is  the  width  of  each  QW,  and  Lq  is  the  thick¬ 
ness  of  the  barrier  between  adjacent  QWs.  Ei  is  the  ground- 
state  energy  of  electrons  at  5^(0  =  0,^  and  is  the 

transmission  coefficient  of  electrons  with  incident  energy  E 
through  a  biased  barrier.^  Because  the  SET  time  is  on  the 
order  of  ^  1  /iS,  which  is  much  smaller  than  the  variation 
time  of  £b(t),  T[E,£i^(t)~\  can  be  evaluated  by  using  the 
time-independent  Schrodinger  equation  for  each  value  of 

4(0. 

When  d£ij(t)ldt¥^0  in  the  nonadiabatic  limit,  there  ex¬ 
ists  a  current  surge  I^(t)  in  addition  to  which  is 

given  by 


S/TTfi^)  lim  (1/Ar)  f  ^  dE{fo[E+Ei 

Af->0  JO 


+  eLB^td£h(t)/dt]-fo[E+Ei-]} 

^-fo[Ey]{m*e^SLB/'rrh^)d£h(l)/dt.  (3) 

When  we  ignore  /^(t),  CDs  in  different  QWs  all  equal  «2i)  • 
The  existence  of  4(0  drives  the  CD  in  each  QW  away  from 
W2Z).  The  charge  fluctuation  A  2(0  itt  the  QW  induces  the 
4a(t)  in  the  system,  which  gives  rise  to  an  additional  nona¬ 
diabatic  current 


^na(0  =  eS{va[£hit)  +  fna(0]neff[4(0  + 

(4) 

In  terms  of  £jJit),  we  can  calculate  AQ{t)  in  the  QW.®  The 
quantum-mechanical  continuity  equation  /na(  0+4(0 
+dAQ(t)/dt=0  leads  to  the  following  dynamical  equation 
for  4a(0: 

^Qw[4a(0]^4a(0/^^ 

-  CQw[0]rf4(0/^t-(^5/L5){i;,[4(0 

+  4a(0>eff[4(0  +  4a(0,4]  ‘ 

-vl£b(t)]nMt)J,]},  (5) 

where  CQw[4a(0]  is  the  QW  capacitance.^^  The  total  tun¬ 
neling  current  is  4[4(0]+4a(0.  If  th®  inequality  |4a(0l 
^|4(0|>A«'c(4)/^I^jB  is  satisfied  in  the  Mmit  of  slowly  vary¬ 
ing  4(^)»  we  can  expand  Eq.  (5)  to  the  first  order  in  4a(0. 
Consequently,  Eq.  (5)  is  simplified  to 

d£n.(t)fdt=^d£,(t)/dt-£,,m^  (6) 

where  Rd[£b(0]  is  the  differential  tunneling  resistance.^^ 
Equation  (6)  is  identical  to  the  previously  derived  circuit 
equation^  by  rewriting  V'^(0=L54(0  nnd  Vc(t)  —  Vi^(t) 
“L54a(0.  Equation  (6)  can  be  viewed  as  the  charging 
/discharging  circuit  equation  with  respect  to  4a(0  in  the 
presence  of  a  source  term  d£ij(t)ldt,  in  which 
i'^^[4(0]^Qw[0]  plays  the  role  of  MCDT.  This  uncovers 
the  microscopic  origin  of  the  circuit  model^  and  explicitly 
relates  the  MCDT  to  the  change  of  microscopic  tunneling 
current  in  the  MQWs. 

III.  NUMERICAL  RESULTS  AND  DISCUSSIONS 

In  our  calculation,  we  assign  two  forms  to  4(0  as  fol¬ 
lows: 

✓ 

2  A5,{l/2+l/irtan-'[(f-r,)/r,]} 

£i>(t)  =  {j=^  (7) 

.  sin[2'n-(f+ro)/7’p]- 

InEq.  (7),  (tj+i  —  tj),  Tj,  and  A5^,  respectively,  denote  the 
delay  time,  transition  time,  and  height  of  the  yth  step  in  the 
step-hke  4(0>  while  4»  4’  27rto/4  represent  the 

amplitude,  period,  and  initial  phase  in  the  sinusoidal  4(0* 
The  sample  we  consider  is  a  MQW  structure  composed  of 
GaAs  wells  and  Alo^GaojAs  barriers.  The  parameters  used 
in  the  calculation  are:  Ly^—5Q  A  ,  Ls—339  A  ,  4=40  K, 
4=2  kV/cm,  i;^=2X10^  cm/s,  n2i>=4.lx  10^^cm~^,  and 
5=2.257rX  10”"^  cm"^.  The  other  parameters  will  be  given 
in  the  text. 

Figure  1  displays  both  the  step-like  4(^)  (dashed  curve) 
and  4(0  +  4a(0  from  Eqs.  (5)  (solid  curve)  and  (6)  (dash- 
dotted  curve)  as  a  function  of  time,  where  the  latter  deter¬ 
mines  4[4(0]  +  4a(0  according  to  Eqs.  (1)  and  (4).  In  the 
calculation,  we  take  and  4= At/20  with  At=50  s 

for  j=  1,2, ...  ,10,  Moreover,  we  choose  A4  =  A4o“  “20 
kV/cm  and  L£j—5  kV/cm  with  7  =  2,3, ...  ,9.  From  Eq.  (6), 
we  find  4a(^)  is  either  positive  or  negative  whenever  4(0 
steps  up  or  down,  due  to  the  4(0  that  results.  This  corre- 
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FIG.  1.  Step-like  5^(0  (dashed  curve)  and  €b(t) 
+  ^iia(0  calculated  from  Eqs.  (5)  (solid  curve)  or  (6) 
(dash-dotted  curve)  as  functions  of  r,  where  the  dotted 
line  corresponds  to  zero  field  and  the  airow  indicates 
the  time  at  which  5^(if)=0.  The  upper  inset  shows 
CQw[^na(0]  (in  units  of  10"^°  F)  as  a  function  of  t  (in 
5),  while  the  lower  inset  presents  log{4[^^(t)] 
+/na(r)} (current  in  A)  as  a  function  of  Si,(t)  (in  kV/ 
cm)  with  arrows  indicating  how  £^^(t)  is  varied 
with  t. 


spends  to  either  discharging  or  charging  of  electrons.  The, 
arrow  indicates  the  time  at  which  5^(r)  =  0.  The  finite 
at  this  time  is  associated  with  the  ZBRTC,^  as  seen  in  the 
lower  inset.  Here,  the  current  ripples  due  to  I^(t)  can  be 
observed  each  time  5^(0  jumps.  The  nonlinear  model  in  Eq. 
(5)  suppresses  the  effect  of  I^it)  by  a  factor  of  up  to 

^  where  is  shown 

in  the  upper  inset.  Cqw[ ^naCO]  develops  a  series  of  peaks  at 
the  time  when  5^(r)  jumps.  The  peak  strength  increases  with 
Rd[^b(^)y  The  current  oscillations  at  ±20  kV/cm  seen  in 
the  lower  inset  of  Fig.  1  ^e  due  to  the  current  bistability.  We 
will  describe  this  bistability  and  the  resulting  current  insta¬ 
bility  when  we  discuss  Fig.  2. 

In  Fig.  2  we  present  the  logarithm  of  /fl[^t(0]+^na(0» 
using  Eqs.  (5)  (solid  curve)  or  (6)  (dashed  curve)  as  a  func¬ 
tion  the  sinusoidal  In  the  calculation,  lO”^  s,  Tp 
=400  s  and  5^;^=  18  kV/cm  are  chosen.  From  the  figure,  we 
see  a  current  arch  supported  by  /a[^^(0]+^na(0==0  at  the 
two  symmetric  finite  values  of  (the  current  “offset”  of 
Ref.  1).  The  nonlinear  model  tends  to  squeeze  the  width  of 
the  current  arch  and  increases  the  ZBRTC.  For  a  sinusoidal 


again  find  the  current  oscillations  around  ±18 
kV/cm  due  to  a  bistability  in  the  current.  The  left  inset  shows 
5^(t)  (dashed  curve)  and  ^^(^)  +  ^na(0  from  Eqs.  (5)  (solid 
curve)  and  (6)  (dash-dotted  curve)  as  a  function  of  time, 
where  the  current  oscillations  are  found  in  around 

d£^{t)ldt=0.  From  it,  we  further  find  the  deformation  of 
f^(r)-h£jia(0  relative  to  £b{t),  which  develops  into  a  phase 
delay  between  them  when  r^=40  s  (not  seen  here).  The 
finite  zero-current  fields  are  indicated  by  the  two  arrows  in 
the  right  inset,  where  the  current  hysteresis  can  be  clearly 
seen. 

The  current  instability  arising  at  high  electric  fields  can 
be  explained  as  follows.  A  current  bistability  is  found  in  this 
system  within  the  nonlinear  model  which  supports  the  solu¬ 
tion  with  d8j^^{t)ldt—0  but  around  d£i,{t)ldt—^. 

Due  to  the  current  hysteresis  as  shown  in  the  right  inset  of 
Fig.  2,  we  can  always  find  a  pair  of  bistable  points  which 
show  the  same  total  tunneling  current  at  and  5^2*  Th^ 
field  difference  -£’£,2!  between  the  bistable  points 

decreases  as  |£i(0l  increases.  From  Eq.  (5)  we  know  when 


FIG.  2.  Results  of  log{4[£’i,(t)]+/„a(0}  using  Eqs.  (5) 
(solid  curve)  or  (6)  (dashed  curve)  as  a  function  of  sinu¬ 
soidal  The  left  inset  shows  Si,(t)  (dashed  curve) 
and  ^^(t)  +  £’jia(f)  calculated  from  Eqs.  (5)  (solid  curve) 
or  (6)  (dash-dotted  curve)  (in  kV/cm)  as  a  function  of  t 
(in  j),  where  the  dotted  line  represents  zero  field.  The 
right  inset  presents  the  calculated  4[^&(0]+/na(0  (in 
units  of  10“^®  A),  as  a  function  of  Ei{t)  (in  kV/cm), 
where  the  two  arrows  indicate  the  finite  zero-current 
fields. 
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Bias  (V) 


HG.  3.  Experimental  data  for  the  logarithm  of  Ihe  tun¬ 
neling  cuirent  as  a  function  of  the  applied  sinusoidal 
bias.  The  left  inset  shows  the  bias  (in  V),  the  curve  with 
the  large  oscillating  amplitude,  and  the  measured  cur¬ 
rent  (in  units  of  10“^®  A),  the  curve  vrith  the  small 
oscillating  amplitude,  as  functions  of  time  (in  s),  where 
the  dash-dotted  line  represents  the  zero  current  and 
bias.  The  right  inset  presents  the  measured  current  (in 
units  of  10“^®  A)  as  a  function  of  the  bias  (in  V),  where 
the  two  arrows  indicate  the  finite  zero-current  biases. 


d£jy{t)ldt’=^0,  dS^{t)ldt  approaches  zero  if  ^^i,=  |£^na(0|  at 
a  particular  time.  This  causes  dS^t)ldt  to  switch  its  sign 
before  and  after  this  time,  leading  to  an  instability  of  the 
tunneling  current.  We  further  note  that  the  current  instability 
can  survive  only  around  kV/cm  where  the  con¬ 
dition  can  be  satisfied  and  will  die  away  as 

S£h>\£,^J^t)\  for  smaller  \£b{,t)\.  However,  in  Eq.  (6)  the 
cuirent  instability  is  absent  and  replaced  by  a  decaying  pro¬ 
cess  for  fnaCO  which  ends  with  ^^na(0/^^=4a(0  =  0. 

IV.  EXPERIMENTAL  RESULTS  AND  COMPARISONS 

The  experiment  was  designed  to  measure  the  tunneling 
current  in  the  MQWs  under  a  sinusoidally  varying  bias  with 
frequency  as  low  as  0.01  Hz.  The  device  was  mounted  in  a 
Dewar  cooled  with  liquid  He,  and  a  darkened  metal  plate 
was  placed  at  the  aperture  of  the  Dewar.  A  cold  shield  was 
inserted  over  the  device.  All  these  precautions  were  taken  to 
ensure  that  we  measured  the  current  in  the  absence  of  pho¬ 
tons.  A  function  generator  was  used  to  generate  the  sinu¬ 
soidal  wave  form,  and  an  electrometer  with  femto-amp  reso¬ 
lution  was  used  to  measure  the  current.  The  electrometer  had 
an  analog  voltage  output  which  was  connected  to  a  digitizing 
oscilloscope  to  display  the  signals  with  frequencies  down  to 
0.01  Hz.  The  applied  and  measured  wave  forms  were  ob¬ 
served  simultaneously  on  the  oscilloscope. 

Figure  3  presents  the  logarithm  of  the  experimentally 
measured  tunneling  current  as  a  function  of  the  applied  bias. 
The  sample  used  in  the  measurement  is  a  50  period 
MQW  structure  and  its  parameters  are 
listed  as  follows:  L^=50  A  ,  1^  =  309  A  ,  r,=40  K,  n2D 
=  2.5X10^^  cm*'^  and  5=10“^^  cm‘^.  We  used  a  sinu¬ 
soidal  bias  with  14  =  6  V  and  7^=40  s,  where  is  related 
to  Vjn  through  5;„=14/50(L5+L^),  The  pillars  of  the  cur¬ 
rent  arch  predicted  in  Fig.  2  are  observed  at  ±5,5  V  in  our 
measurement.  The  left  inset  proves  the  existence  of  the  cur¬ 
rent  instability,  shown  as  an  additional  peak  superposed  on 


the  shoulder  of  the  sinusoid-like  tunneling  current  with  a 
phase  delay  with  respect  to  the  bias.  The  right  inset  confirms 
the  current  hysteresis  predicted  in  Fig.  2.  The  fattened  shape 
of  the  current  hysteresis  curve  makes  it  very  difficult  to  ob¬ 
serve  the  current  oscillations  with  more  than  one  period  since 
the  field  difference  between  the  bistable  points  is  very  large. 
The  current  ripples  seen  in  Fig.  1  are  also  observed  from  our 
experiment  and  will  be  presented  elsewhere. 

V.  CONCLUSIONS 

In  conclusion,  by  proposing  the  CSM,  we  have  explored 
the  mystery  of  the  observation  of  the  microscopic  tunnehng 
current  change  on  a  macroscopic  time  scale  and  found  the 
microscopic  origin  of  the  previously  reported  circuit  model 
used  to  explain  the  ZBRTC.  The  physics  of  the  tunneling 
process  with  a  very  large  time  constant  in  the  MQWs  has 
been  elucidated.  With  the  CSM  we  have  predicted  some  ad¬ 
ditional  features  of  the  tunneling  current  in  the  presence  of  a 
time-dependent  bias,  including  a  current  ripple,  current  hys¬ 
teresis,  a  current  arch,  and  current  instability.  All  of  these 
predictions  have  been  confirmed  by  our  experiment. 
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